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Abstract 

The dynamics of two nonlinear Bloch systems is studied from the viewpoint of bifur- 
cation and a particular parameter space has been explored for the stability analysis 
based on stability criterion. This enables the choice of the desired unstable periodic 
orbit from the numerous unstable ones present within the attractor through the pro- 
cess of closed return pairs. A generalized active control method have been discussed 
for two Bloch systems arising from different initial conditions. 
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1 Introduction 



The simple Bloch equations exhibit the dynamics of an ensemble of spins which 
do not exhibit mutual coupling. The magnetization dynamics can be determined 
from the linear nature of these equations. However, in high field Nuclear Magnetic 
Resonance(NMR), these simple structures break down due to the presence of radiation 
damping[l] or a demagnetizing field [2]. Manifestations of nonlinear spin dynamics 
could be observed due to the presence of an additional field which is proportional to 
the components of the magnetization. The dipolar magnetizing field could be shown 
to give rise to multiple echoes in liquid helium [2] and in water at high magnetic 
field [3]. Moreover, radiation damping or demagnetizing field effects [4-6] can result 
in the existence of pseudo-multiquanta peaks. Complex nonlinear behaviour of the 
nuclear spins at high magentic field was recently demonstrated through manipulation 
of a large magnetization by radiation damping based electronic feedback and giving 
rise to a train of steady state maser pulses of the water magnetization [7]. 

Recently, Abergel[8] investigated the possibility of observing chaotic solutions 
of the Bloch equations. The various anomalies that arises in NMR experiments have 
been studied in terms of chaos theory. Ugar et al extended the calculation of Abergel 
and demonstrated the synchronization of two such Bloch systems through 'Master' 
and 'Slave' arrangement through a suitably designed controller. 

The work of OGY and Pecora and Carroll [9] led to wide applications outside 
the traditional scope of chaos and nonlinear dynamics. This has led to the estab- 
lishment of two active areas of research viz. synchronization and control. Recent 
problems casts the problem of chaos synchronization in the framework of nonlinear 
control theory. This unifies the study of chaos control and chaos synchronization. 

Pyragas have proposed two methods of permanent chaos control with small 
time continuous perturbation in the form of linear feedback[10]. The stabilization of 
a unstable periodic orbit (UPO) of a chaotic system is achieved either by combined 
linear feedback with the use of a specially designed external oscillator or by delayed 
self controlling linear feedback. 

An open-plus-closed-loop (OPCL) method of controlling nonlinear dynamic 
systems was presented by Jackson and Grosu[ll]. The input signal of their method 
is the sum of Hubler's open-loop control and a particular form of a linear closed- 
loop control, the goal of which can be selected as one of the UPO's embedded in the 
chaotic attractor or another possible smooth functions of l t\ The asymptotic stability 
of the controlled nonlinear system is realized by the linear approximation around the 
stabilized orbit. But the calculation of the closed-loop control signal is very difficult 
in some cases, specially for complex and high dimensional chaotic systems. 

Yu. et al [12] proposed a method for controlling chaos in the form of special 
nonlinear feedback. The validity of this method based on the stability criterion of 
linear system and can be called stability criterion method(SC method). The con- 
struction of a nonlinear form of limit continuous perturbation feedback by a suitable 
separation of the system in the SC method does not change the form of the desired 
UPO. The closed return pair technique[13] is utilised to estimate the desired periodic 
orbit chosen from numerous UPO's embedded within a chaotic attractor. The advan- 
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tage of this method is that the effect of the control can be generalized directly without 
calculation of the maximal lyapunov exponent of the UPO using the linearization of 
the system. This method does not require linearization of the system around the 
stability orbit and calculation of the deviation at UPO's. The method has been used 
by researchers in the control of Rossler system, chaotic altitude motion of a spacecraft 
and the control of two coupled Duffing oscillators. 

In this communication, the bifurcation analysis of the nonlinear Bloch equa- 
tions in a particular parameter regime has been discussed and stability criterion 
method has been used to investigate the synchronization of two such Bloch systems. 
The technique of generalized active control is also used for deriving a controlled tra- 
jectory of the system. 



2 Formulation 

The model is derived from a magnetization M precessing in the magentic induction 
field B in the presence of a constant radiofrequency field B x with intensity Bi = ^ 
and frequency u> r f. A magnetization-dependent field B fB [14]: 

B FB = 7GM t e-^, 

with 

M t = M x + iM y 

where G is the enhancement factor with respect to the magnitude of the magnitude 
of the transverse magnetization and ip is the phase of the feedback field. The following 
modified nonlinear Bloch equation govern the evolution of the magnetization, 

M x = 5M y + GM Z (M X sin(^) - M y cos(^)) - ^ (2.1) 

J-2 

M y = -5M X - u x M z + GM Z (M X cos(^) + M y sin(^)) - ^ (2.2) 

-*2 



M z = u\M y - Gsim»(il4 2 + My) - (Mz Ti M ° ; (2.3) 



where 5 = u r f — u and T 1; T 2 are the longitudinal time and transverse relaxation 
time respectively. The above three equation is transformed by introducing the reduced 
dimensionless variables: 

G 

t — > ujit, G — > M — = 7, 5 — > 5/cui, Ti 2 — > oj\Ti 2 and 

M, -> M x /M ./• M, :l - M y /M n M, — M z /M = z. 
The nonlinear Bloch system is then expressed as, 

x 

x = 5y + jz(x sin(c) — y cos(c)) — — (2.4) 

I 2 
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y = —5x — z + r yz(x cos(c) + y sin(c)) — — (2.5) 

z = y- 1 sm(c)(x 2 + y 2 )-^— ^ (2.6) 

1 1 

2.1 Analysis of the chaotic Dynamics 

The presence of the attractor has already been shown in the previous analysis of the 
system. The system possess two attractors for two different set of parameters viz., 

7 = 35.0, 5 = -1.26, c = 0.173, r\ = 5.0, T 2 = 2.5 

and 

7 = 10.0, 5 = 1.26, c = 0.7764, T t = 0.5, T 2 = 0.25 

The behaviour of the system with a change in the parameter 7 is shown with the 
help of bifurcation diagram. We have shown the variation of the parameter from 
7 = 23.0 to 7 = 32.0 . The dynamics of the system is shown in Figure(l) where the 
system exhibits multiperiodicity for higher values of the parameter 7 whereas it is 
clear that the chaoticity arises from 7 = 30.0 and extends upto 7 = 37.0. In the range 
of 32.5 < 7 < 33.5, there is a sudden transition to periodicity. The unstable periodic 
orbit finally attains chaos which continues upto 7 = 37.0 and after that again the 
system intermittently transits to multiperiodic regime. 

A detailed stability analysis based on largest lyapunov exponent of the sys- 
tem was carried out over a certain parameter region. The two parameters 7 and c 
were taken into consideration and the system dynamics was distinguished into three 
different categories according to value of the largest positive lyapunov exponent. 

In the stability plot the range of the x-axis and the y-axis are 22 < 7 < 32and 
—0.2 < c < 0.7 respectively. The regions representing the dynamics of the system in 
the (7, c) plane is shown in Figure(2). The regions marked by cross ' x ' depicts chaotic 
state where largest positive lyapunov exponent A > 0.05. A dark dot '•' indicates 
0.003 < A < 0.05. These regions exhibit multiperiodic as well as chaotic behaviour. 
Finally, the rest of the plane which is covered by plus sign '+ 'corresponds to purely 
stable behaviour with A < 0.003. The system illustrates a periodic behaviour or 
equilibrium state under these parameter conditions. 



2.2 Control based on stability criterion 

To control the system of nonlinear Bloch equation a time continuous nonlinear dy- 
namic system with input perturbation described by 

^ = /(*(*))+«(*) (2.7) 

where x G R n and u G R n are the state vector and input perturbation of the system 
respectively. Equation (2.7) without input signal (u = 0) has a chaotic attractor Q. 
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A mapping / : R n — > R n is denned in n- dimensional space. We suitably decompose 
the function f(x(t)) as, 

f(x(t))=g(x(t)) + h(x(t)) (2.8) 

where the function g(x(t)) = Ax(t) is suitably disposed as a linear part of f(x(t),t) 
and it is required that is a full rank constant matrix, all eigenvalues of which ahve 
negative real parts. So the function h(x(t)) = f(x(t)) — Ax is a nonlinear part of 
f(x(t)). Then the system (2.7) can be rewritten as 



dx 

— = Ax(t) + h(x(t))+u(t) 

(Jjv 



(2.9) 



Let D(x(t)) = —h(x(t)), then the function f + D = f — h, is a linear mapping 
with respect to the state vector x, namely, 



(f + D)(x)=Ax 



(2.10) 



Let, x*(t) = x*(t + jT), j = 1, 2, ... be a period- j trajectory embedded within 
Q. The input signal u(t) is considered as a control perturbation signal as follows 

u(t) — D(x(t)) — D(x*(t)) (2.11) 

Substituting Eq.(2.11) into Eq.(2.9) system (2.7) and (2.9) can be rewritten as, 

x-x* = (f + D)(x) - (/ + D)(x*) = A(x - x*) (2.12) 

The difference between x(t) and x*(t) is defined as an error w(t) = x(t)—x*(t), 
then evolution of which is determined by Eq.(2.12) as, 



w(t) = Aw(t) 



(2.13) 



Obviously, the zero point of w(t) is its equilibrium point. Since all eigenvalues 
of the matrix A have negative real parts, according to the stability criterion of linear 
system, the zero point of the error w(t) is asymptotically stable and w(t) tends to 
zero when t — > oo. Then the state vector x(t) tends to the period-j trajectory x*(t). 
It implies that the unstable periodic orbit is stabilized. Stable soltions belonging 
to different basins of initial conditions can also be the alternative solutions of very 
complicated periodically driven dynamic systems along with the stabilized UPO. The 
stabilization is obtained by modifying Eq.(2.11) as follows 

u{t) = D(x(t)) - D(x*(t)) 

= A(x - x*) + f(x*) - f(x), if | x - x* \< e 

= 0, otherwise 

For our system of equations given by (2.4-2.6) we can write the linear and nonlinear 
part as, 



/(*) 



-5 



V o 



r 2 



\ 



-1 



( x\ 



V z J 



( 7^(x sin(c) — y cos(c)) ^ 

•yz(x cos(c) + y sin(c)) 
^ -7sin(c)(a; 2 + y 2 ) + ± J 



(2.14) 
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where 



A = 



"A 6 N 



-i 



V o 



h(x) = 



( ^yz(x sin(c) — y cos(c)) \ 

•yz(x cos(c) + y sin(c)) 
^ -7sin(c)(a; 2 + y 2 ) + £ J 



Then the control perturbation signal is derived as follows 

/ — 7^(xsin(c) — y cos(c)) + r )z*(x* sin(c) — y* cos(c)) \ 



u{t) = 



— ^yz(x cos(c) +ysin(c)) + jz* (x* cos(c) + y* sin(c)) 
y 7sin(c)(x 2 + y 2 ) — 7sin(c)(a;* 2 + y* 2 ) j 



In order to obtain the necessary information on an appropriate location of a desired 
periodic solution x*, the strategy of the close return pair is taken into account. The 
method is to generate a time series of the chaotic trajectory by stroboscopically 
sampling in every period T when u = 0. 

The data sampling can be used to detect the close return pairs, which consists of two 
successive points nearing each other, and indicate the existence of a periodic orbit 
nearby. Let x it i and x i>2 are used to denote thr first point and its successive point 
of the ith collected return pair, = 1,....,M respectively, where M is the maximum 
number of collected return pairs. When the first close return pair has been detected, 
taking the first point reference point, a number of close return pairs nearing 

the reference point can be detected. 



x i,l — x l,l |^ Cl ? | x i,2 — x l,l |^ e 2) 



1,2, 



M 



We define the mean value as, 



1 



M 



x ' = 2M&i 1 i + I <>2] 

where x* is regarded as an appropriate fixed point. This fixed point can be used to 
define a restriction condition | x(t) — x*(t) \< e within which the control input signal 
u ^ 

We have targeted both the attractors as described in the previous section. In 
both cases M = 3. In the first case, x*(t) = -0.221, y*(t) = -0.021, z*(t) = 
0.141 and e = 0.45 and in the second case, x*(t) = 0.219, y*(t) = -0.316, z*(t) = 
0.790 and e = 0.40. Figures (3a) and (3b) respectively show the results of control. 
After a short transition both the attractor stabilizes on a periodic trajectory. 



2.3 Synchronization using generalized active control 

Bai and Lonngren[15] proposed an active control process and two Lorenz systems were 
synchronized using that technnique. Ho and Hung used the method of generalized 
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active control to synchronize two completely different systems[16]. Here we show the 
synchronization of two nonlinear Bloch equations and this technique is different from 
the one by Ugar et al[17] in the sense that the control signals in this particular case 
do not contain any positive gain term. It is found to be equally effective and the 
results can be compared to that obtained earlier. A second nonlinear Bloch equation 
is considered with the same parameter values as the previous one but differing in the 
initial conditions. 

u 

ii = 8v + 7w(m sin(c) — v cos(c)) — — + r]i(t) (2-15) 

r 2 

V 

v = —5u — w + ^w{u cos(c) + v sin(c)) — — + 772(0 (2-16) 

I 2 

w = v - 7 sin(c)(w 2 + v 2 ) - ^~ - + Vs(t) (2.17) 

1 1 

Subtracting Equations(2.15-2.17) from (2.4-2.6) and performing the required calcu- 
lations, the control signals are obtained as follows 

r)i(t) = (i(t) — ^w{u sin(c) — v cos(c)) + jz(x sin(c) — y cos(c)) (2-18) 

V2{t) — C2W — 7w(mcos(c) + i>sm(c)) + 7,2(2; cos(c) +ysin(c)) (2-19) 

m(t) = (s(t) + 7 sin(c)(M 2 + v 2 ) - 7sin(c)(a; 2 + y 2 ) (2.20) 



where 



1 £1 - 5e 2 



£ 2 (t) = 8e 1 + (— - l) e 2 + e 3 



Cs(0 = - l) £3 - £2 

with E\ = u — x, e 2 = v — y, 63 = w — z. The parameters are 7 = 35.0, 5 = 
— 1.26, c = 0.173, T 1 = 5.0, T 2 = 2.5 and the error dynamics is shown in Figure (4) 
which converges zero indicating the synchronization between the two systems which 
evolve from two different states. 



3 Conclusion 

In our above analysis we have studied a different methodology of synchronization in 
a Maxwell-Bloch system by analyzing its bifurcation pattern and stability. The main 
emphasis is on identifying a unstable periodic orbit and to adopt the strategy of closed 
return pairs to control whose effectiveness is represented in our results. Regarding 
the method of active control, while the previous authors were concerned about two 
systems with different parameter values our consideration is on two different initial 
state of the system-that is their magnetization. Our approach does not take into 
account any gain term as it is more generalized than the previously considered one. 

Acknowledgement:: One of the authors B.Rakshit is thankful to CSIR, Govt, of 
India for a research fellowship in a project. 
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5 Caption for Figures 

[1]. The bifurcation diagram of the Bloch system with respect to the parameter 7. 
[2]. Chaotic domain in the (7,0) plane. '•' and '+' represent stable states whereas 
'x' denotes chaos. 

[3]. Temporal evolution of the system after the application of the control perturba- 
tion due to stability criterion, a) 7 = 35.0,5 = — 1.26, c = 0.173, T\ = 5.0, r 2 = 2.5; 
b)7 = 10.0, 5 = 1.26, c = 0.7764, n = 0.5, r 2 = 0.25; . 

[4]. Error dynamics of the system due to control by generalized active control. 
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